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' In this paper, we have reviewed the present status of the theory of equilibrium config- 

, urations of compact binary star systems in Newtonian gravity. Evolutionary processes of 

jy) ■ compact binary star systems due to gravitational wave emission can be divided into three 

stages according to the time scales and configurations. The evolution is quasi-stationary 
until a merging process starts, since the time scale of the orbital change due to gravitational 
wave emission is longer than the orbital period. In this stage, equilibrium sequences can 
be applied to evolution of compact binary star systems. Along the equilibrium sequences, 
there appear several critical states where some instability sets in or configuration changes 
drastically. We have discussed relations among these critical points and have stressed the 
importance of the mass overflow as well as the dynamical instability of orbital motions. 

Concerning the equilibrium sequences of binary star systems, we have summarized classi- 
cal results of incompressible ellipsoidal configurations. Recent results of compressible binary 
star systems obtained by the ellipsoidal approximation and by numerical computations have 
been shown and discussed. It is important to note that numerical computational solutions 
to exact equations show that compressibility may lead realistic neutron star binary systems 
to mass overflows instead of dynamical disruptions for a wide range of parameters. 

o 

§1. Quasi-equilibrium stage of evolution of compact binary star systems 

1.1. Time scales and classification of evolutionary stages of compact binary star 
systems 
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As is well known, gravitational waves from compact binary star systems will 
carry away the angular momentum and energy from the systems. It implies that 
the binary star system will necessarily evolve to shrink its orbital radius and that 
eventually merging of two compact component stars will occur. There are three 
typical time scales to classify such an evolution: 1) the rotational period, r rot , 2) the 
time scale of the orbital change due to gravitational wave emission, tgw> an d 3) the 
light crossing time scale of the system, Td yn - 

Let us consider a binary system consisting of two stars with masses M\ and M<i . 
Two stars are in a circular orbit with a separation of A = 2d. Here the separation 
of the binary system is defined by the distance between two centers of mass of two 
stars. If stars are point masses with the same mass and the quadrupole formula 
for the gravitational radiation reaction is used, time scales mentioned above can be 
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estimated as follows in Newtonian gravity: 



Tdyn = d , (1-1) 

/ 1 dA\- x 5 /Mr 3 / d \ 3 
^W^hr-s- «o -5" "5- > C 1 - 2 ) 



A dt J 8 V#*/ V# 

^T^UJ UJ ' (1 " 3) 

where i?* and 17 are the radius of the star and the angular velocity, respectively, and 
geometrized units c = G = 1 are used. If we specify the value of the compactness 
M/R*, we can evaluate these time scales as shown in Table I and Fig. 1. In Table I 
ratios of time scales, the rotational velocity v mt and the number of rotations during 
the time scale tgw> N, are shown for two cases of the compactness M/R*. In Fig. 1, 
the ratios of time scales are plotted against the separation normalized by the stellar 
radius. 

Table I. Ratios of time scales, the rotational velocity v vo t and the number of rotations, N, during 
tgw at the distance of several stellar radii for different compactness parameters. 
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Fig. 1. Ratios of time scales logfYow/Tdyn) (thick curves) and log(T rot /T d yn) (thin curves) are plot- 
ted against the separation normalized by the stellar radius for two values of M/R,. M/R, =0.1 
(solid curves) and 0.2 (dashed curves). 



As seen from Table I and Fig. 1, since the time scales vary considerably as a 
function of the separation, we can divide the evolution of binary star systems due to 
gravitational wave emission into three distinct stages: 1) the quasi-equilibrium stage 
where tgw ^ Trot or d ^ i?*, 2) the merging stage where tqw ~ Trot or d ~ i?*, 
and 3) a neutron star or a black hole formation stage where at first row ~ Trot or 
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d ~ R* and later tgw ^ Trot or d & i?*. In stage 1, gravitational wave emission 
hardly affects the orbital motion during several rotational periods, while in stages 
2 and 3 an appreciable change of the orbital motion and/or configuration occurs in 
one or less orbital period or on a dynamical time scale. 

Therefore, hydrodynamical simulations are most appropriate for investigations of 
stages 2 and 3 but not for stage 1. For stage 1, it is proper to assume that binary stars 
are in quasi-equilibrium states because there is little change of the configuration and 
the orbital motion of the binary star system. Moreover, we need to take the structure 
of the stars into account for the state with rather small separations d/R* ^2 where 
this quasi-equilibrium approach may be still applied. Evolution of a binary star 
system in this stage is well approximated by an equilibrium sequence of a binary 
with a finite size, if we take account of proper constraints in construction of the 
sequence. 

1.2. Evolution of velocity fields and equilibrium sequences 

One equilibrium state of a binary star system can be specified by the following 
quantities: 1) masses of the stars, 2) the value of the separation, 3) the equation of 
state and 4) rotation laws of the stars including the internal motions. 

Concerning the first two, we can freely choose those values. Although the equa- 
tion of state for the realistic neutron star matter has not been fully understood for 
the high density region above the nuclear density, we will be able to try several sets 
of equations of state. 

However, it is not easy to specify the rotation law for binary star systems because 
the rotation law depends on the initial state, the physical state of the matter, the 
evolutionary time and so on. In this paper, we will consider two extreme cases for 
viscosity: 1) strong viscosity and 2) weak viscosity. If the effect of viscosity is strong 
enough, the internal flow would synchronize with the orbital rotation. On the other 
hand, if it is weak, the initial state of the flow would be maintained. 

As KochanekEP and Bildsten and Cutler a have shown, the situation for neutron 
stars becomes very simple. The effect of viscosity for realistic neutron stars is so small 
that it may not affect the rotation law of binary neutron star systems. This implies 
that the vorticity of the neutron star system is conserved during evolution because 
the radiation reaction force due to gravitational wave emission is the potential force 
at the leading order. 

Furthermore, since values of the angular velocity and the vorticity, £, of neutron 
stars will increase as the orbit shrinks, the following relations will be established at 
the final state of the quasi-equilibrium stage: 



where quantities with suffix j or t denote those at the initial state or at the final 
state, respectively. Here Q and are the vorticities in the rotating frame. Since the 
vorticity in the inertial frame Co is conserved, we can obtain the following relation: 



Q f \ > \Oi\, 
IC/I » 




(1-5) 



ICol « 1%!, IC/I 
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This means that the vorticity in the inertial frame can be considered to be neg- 
ligibly small compared with the angular velocity and the vorticity in the rotating 
frame at the final state of the quasi-equilibrium evolution stage. In other words, we 
may consider that flow fields of binary neutron star systems can be regarded to be 
irrotational. 

Therefore, in this paper, we will treat two extreme cases: 1) synchronously 
rotating binary star systems and 2) irrotational binary star systems. 

In principle, we can obtain equilibrium sequences by connecting equilibrium 
states with the specified rotation law. However, in order to follow quasi-stationary 
evolution of binary star systems due to gravitational wave emission by making use 
of a quasi-equilibrium sequence, we need to consider conserved quantities during 
evolution. For binary neutron star systems, one conserved quantity is of course the 
vorticity in the inertial frame as mentioned above. For binary star systems emitting 
gravitational waves, the baryon number or the baryon mass is the second conserved 
quantity because matter does neither flow out from nor flow into the system. There- 
fore, we can construct equilibrium sequences by connecting equilibrium states with 
the same baryon mass and use them as quasi-stationary evolution sequences due to 
gravitational wave emission. 

1.3. Characteristic feature of equilibrium sequences 

Equilibrium sequences can be characterized by several sets of physical quantities 
such as the total angular momentum - separation relation, the angular velocity - 
separation relation, the energy - separation relation and so on. In this paper, we 
will choose the total angular momentum - separation relation as the typical relation 
which reflects the essential feature of the equilibrium sequence or, in other words, of 
the quasi-stationary evolution sequence of the binary star system due to gravitational 
wave emission. 

In order to show the characteristic feature of the equilibrium sequence, we will 
assume that masses of two stars are the same and that the stars rotate in a circular 
orbit. The total angular momentum J tot is divided into two parts as follows: 



where </ rb and J sp i n are the orbital angular momentum and the spin angular mo- 
mentum, respectively. If we assume that the stars are congruent, they are expressed 
as: 




spin j 



(1-7) 
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where we have used the relation 



n 2 A 3 = 2M 
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Here /, p and v v are the moment of inertia, the density and the (^-component of 
the internal velocity of the component star in the rotating frame, respectively, and 
(r,9,ip) are the coordinates whose origin coincides with the center of mass of the 
component star. 

From these expressions, we can see the dependency of the angular momentum on 
the separation as well as on the internal structure of the star. First, as the separation 
decreases, the orbital angular momentum decreases as A 1 / 2 , while the spin angular 
momentum increases as A~ 3 ^ 2 . Second, the internal structure of the finite-size star 
is reflected in the second and third terms of Jtot through / and v^. The moment of 
inertia depends on the matter distribution which is governed by its self-gravity, the 
tidal force from the companion star, the orbital motion, the internal flow including 
rotation and the equation of state. The internal velocity v v is determined from the 
initial state and evolution of the velocity field. 

The first feature implies that there can exist a minimum angular momentum 
state for a certain value of the separation, say A ni { n . From the second feature the 
effect of the internal structure changes the spin angular momentum but keeps the 
orbital angular momentum unchanged. In particular, the change of the moment of 
inertia will result in the shift of the value of A m i n because the moment of inertia 
affects through the term which behaves as A~ 3 ^ 2 . Thus, if / increases (decreases), 
A m i n becomes larger (smaller). The effect of the internal flow is not so drastic. If we 
consider the velocity field which flows in the counter direction with respect to the 
orbital motion, i.e. v v < 0, the spin angular momentum is decreased almost evenly 
for all separations so that the total angular momentum decreases but the position 
of A m \ n is hardly shifted. These features are shown in Figs. 2 and 3. 
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Fig. 2. Schematic figure of the angular momentum J - separation A relation. The effect of the 
moment of inertia is shown. Thin solid curves correspond to J or b and J sp m and thick solid curve 
to Jtot for a certain value of /. Dashed curves correspond to those for larger I. 



§2. Critical states along the equilibrium sequence of binary star systems 

When we consider evolutionary sequences of binary star systems, it is important 
to know whether each model along the sequences is stable or not against a certain 





kind of perturbations. If some kind of instability sets in at some point, equilibrium 
configurations beyond that point do not exist in realistic situations. Even if insta- 
bility does not occur, at some stage there may occur some kind of significant change 
of configurations. In this section we will discuss such kind of critical states which 
will appear along evolutionary sequences. 

2.1. Contact state for congruent stars 

For congruent stars, if the separation decreases, it is most probable that two 
binary stars come to contact by considering the symmetry of configurations in equi- 
librium states. In order to find this critical state, we need to treat deformation of 
configurations very accurately. We will use A c or d c for the separation of this con- 
tact phase or the half of the separation, respectively. After this contact phase, its 
configuration is no more a binary state but a single body is formed. 

2.2. Roche lobe filling state (cusp formation state) 

If masses of two stars are different or if configurations of two stars are different 
such as a black hole - neutron star system or if there exist internal motions, decrease 
of the separation results in a Roche lobe filling state. Here we use the word "Roche 
lobe" in its generalized sense because models can include internal flows. Since there 
is a cusp point on the surface of the Roche lobe, a cusp shape appears on the stellar 
surface, too. This state is denoted by the separation ^of or the half of the separation 
doF- Here the subscript of denotes that the mass begins to overflow from the Roche 
lobe filling star to the other star. 

2.3. Roche limit state (minimum separation state) 

The Roche limit is defined as the minimum separation for a binary system. 
Within the Roche limit there are no equilibrium states for the component stars. 
For a spherical star, the Roche limit state coincides with the state at the maximum 
angular velocity. However, if we take deformation of stars into account, this is no 
more the case. For some binary configurations, equilibrium states with larger angular 
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velocity than that of the Roche limit state can exist. The Roche limit and its half 
distance are denoted as ^4rl and dm,, respectively. 

2.4. Turning point of equilibrium sequence 

As discussed before, along the equilibrium sequence of binary star systems, there 
may appear a state where the total angular momentum becomes minimum. Since 
gravitational waves carry away the angular momentum from the system, the min- 
imum state of the angular momentum denotes a critical state for equilibrium se- 
quences. For synchronously rotating binary star sequences, this state corresponds 
to the onset of secular instability due to viscosity. On the other hand, for irrota- 
tional binary star sequences, dynamical instability of the orbital motion sets in at 
this point. We will use A^ yQ and d^ JQ for the separation and the half of it at this 
state, respectively. 

2.5. Dynamical instability against collapse 

For compact stars such as white dwarfs or neutron stars, we need to consider 
the stability of a single star to collapse. This instability sets in at the state of the 
maximum gravitational mass along the equilibrium sequence with the same distribu- 
tion of the specific angular momentum for barotropes in Newtonian gravity. @) (In 
general relativity situation is somewhat different. &) The criterion of this instability 
can be written as follows: 



where p mauX and /9 m ax,c are the maximum density of the star and the density for the 
maximum gravitational mass state, respectively. It is not easy to find this critical 
configuration exactly because the equilibrium sequences for compact binary star 
systems are not those with the distribution of the specific angular momentum kept 
constant. However, roughly speaking, if the maximum density of the star increases 
as the separation decreases, there arises a possibility of onset of this instability for 
each component star. 

2.6. General relativistic instability 

In this paper, although we treat equilibrium configurations only in Newtonian 
gravity, we should point out that there is a critical state due to general relativistic 
instability of the orbital motion, i.e. the innermost stable circular orbit. The critical 
separation or the distance for this instability is denoted as ^4gr ° r ^GR- Here 



where M% t is the total mass and A is a positive quantity which depends on the 
general relativistic structure of the system. 

2.7. Classification of binary equilibrium sequences by critical states 

By considering the order of appearance of these critical points along the equilib- 
rium sequence, we can classify the quasi-stationary evolution of binary star systems 
into three types with several subtypes from the topology of the equilibrium sequence 
on the total angular momentum - separation plane. 



Pmax ^ Pmax,c 




d GR = 6M tot + A 
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Type la: The equilibrium sequence terminates at either A c or Aop. There is no 
turning point. 

Type lb: The equilibrium sequence terminates at Arl. There is no turning point. 

Type 2: The equilibrium sequence terminates at either A c , Aop or Arl but the 
first critical point along the equilibrium sequence is the turning point. Thus the 
following condition is satisfied: Ad yn > A c , Aop, A^. 

Type 3: The equilibrium sequence terminates at either A c or Aof- The first critical 
point along the equilibrium sequence is the turning point. Thus the following 
condition is satisfied: ^4dyn > ^4rl- 

For quasi-stationary evolution of Type la, the final outcome is a single star or 
mass overflow from one component star to the other. For stars of Type lb, 2 and 
3, they suffer from dynamical disruption eventually. Of course, if the turning points 
are those for secular instability, stars will excite internal motion due to viscosity and 
can come closer. However, at the end they will finally encounter the turning point 
where dynamical instability sets in and suffer from dynamical disruption. Schematic 
figures of these types are shown in Fig. 4. Almost the same analysis has been done 
by Lai, Rasio and Shapiro. U 




A A A 

Fig. 4. The total angular momentum is plotted against the separation for equilibrium sequences 
which belong to Type 1 (left panel), Type 2 (middle panel) and Type 3 (right panel). The Type 
1 equilibrium sequence terminates at either the contact state, the Roche lobe overflow state 
(Type la) or the Roche limit state (Type lb). The Type 2 equilibrium sequence terminates at 
the contact state, the Roche lobe overflow state or the Roche limit state but there is a turning 
point along the sequence. The Type 3 equilibrium sequence terminates at the contact state or 
the Roche lobe overflow state. There is a turning point along the sequence. The Roche limit 
state comes between these critical configurations. 



§3. Classical results for binary star systems 

Until early of 1980s, it was hard to obtain binary configurations of compress- 
ible gases. Therefore, structures of binary star systems had been investigated by 
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assuming that 1) the density is constant, that 2) the shape of the star is ellipsoid, 
that 3) the angular velocity is constant and that 4) the vorticity is also constant. 
Under these assumptions three kinds of classical ellipsoids for binary star systems 
have been solved (see e.g. Ref. ||)). 

3.1. Roche ellipsoid: synchronously rotating sphere (or point mass) - ellipsoid sys- 
tem 

Let us assume that one of the stars is a point mass or a sphere and that the 
other star is synchronously rotating with the orbital motion. If we include only 
the leading term of the tidal force from the point mass or the spherical star, the 
equilibrium condition for the fluid star can be expressed as 

P - = -01 + ^ (* 2 - \y 2 - \z 2 ) + \n\x> + y") + constant , (3-1) 

where p, <pi, are the pressure, the gravitational potential of the fluid star (self- 
gravity) and the mass of the spherical star, respectively. The Cartesian coordinates 
(x, y, z) are used. If the shape of the fluid star is ellipsoid, this equilibrium condition 
can be satisfied by choosing the angular velocity Q appropriately. This ellipsoid is 
called the Roche ellipsoid. 

For the Roche ellipsoids with M\ = M2, the angular momentum - separation 
relation is shown in Fig. 5. As seen from this figure, the Roche ellipsoidal sequence 

3 1 . ■ . ■ . 1 
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Fig. 5. Angular momentum is plotted against the separation for Roche ellipsoids. The orbital 
angular momentum (dashed curve), the spin angular momentum (dash-dotted curve) and the 
total angular momentum (solid curve) are shown. The turnaround of the spin curve is due to 
changes of the internal structure. 

belongs to Type 3. The curve of the spin angular momentum shows a turnaround due 
to a significant change of the internal structure. Along this sequence there appears 
a turning point where secular instability sets in. At this turning point the shape of 
ellipsoid is specified by 02/01 = 0.65 and 03/01 = 0.59 where (01,02,03) are three 
principal axis distances of the ellipsoid. 
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Fig. 6. Same as Fig. 5 but for Darwin ellipsoids. There is no turnaround of the spin curve. 

3.2. Darwin ellipsoid: synchronously rotating congruent ellipsoids 

If two fluid stars are congruent and synchronously rotating, the equilibrium 
condition can be expressed as follows: 

- = -0i + {(5 l x 2 + [3 2 y 2 + l3-iZ 2 ) + \o 2 {x 2 + y 2 ) + constant , (3-2) 
P ^ 

where f3's are constants which can be determined from the leading term of the tidal 
force from the secondary star. For this case, the ellipsoidal configuration satisfies the 
above equilibrium condition. Solutions to this equation are called Darwin ellipsoids. 
The equilibrium sequence is shown in Fig. ^. The sequence terminates at the contact 
state. Thus the Darwin sequence belongs to Type 3 sequence. The Roche limit state 
is characterized by 02/01 = 0.62 and 03/01 = 0.56. 

3.3. Roche- Riemann ellipsoid: a sphere (or point mass) - ellipsoid system with 
vorticity 

Fluid binary stars with nonzero vorticity were first investigated by Aizenman.H* 
The system is almost the same as the Roche problem but in the fluid star there exist 
internal motions represented by the vorticity. The equilibrium condition is written 
as follows: 

- = -0i - 2 + \tt 2 (x 2 + y 2 ) + \u 2 + (2f2 + ()& + constant , (3-3) 

where 02, u and & are the gravitational potential from the spherical star, the velocity 
of the internal motion in the rotating frame and the stream function of the internal 
motion, respectively. This condition can be satisfied by ellipsoids and ellipsoidal 
solutions are called Roche- Riemann ellipsoids. According to the results of Aizenman, 
by introducing the internal flows, the minimum separation, i.e. Roche limit, can 
become smaller than hat for the Roche ellipsoids. In other words, two stars can 
come closer each other. 
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§4. Recent results for binary star systems in Newtonian gravity 

Although basic properties of binary star systems have been understood by us- 
ing classical results mentioned in the previous section, several crucial things have 
remained to be investigated. First, the constant density assumption is not valid 
for realistic binary stars. Second, the real shape of deformed stars is not ellipsoidal. 
Therefore we need to take these points into account to understand realistic evolution 
of compact binary star systems. 

4.1. Ellipsoidal approach 

The effect of the compressibility is introduced by Lai, Rasio and Shapiro. B They 
have investigated polytropic binary star systems: 

p = Kp 1+1/n , (4-1) 

where K and n are a polytropic constant and the polytropic index, respectively. 
Concerning the shape of binary stars, they have assumed that it is approximated 
by ellipsoidal shape. By using the energy variational principle, they have obtained 
approximated equilibrium configurations. 

Their main results are as follows. For simplicity we will call Roche or Roche- 
Riemann ellipsoids for the systems consisting of a spherical star (or point mass) and 
a compressible ellipsoid without or with internal motions, respectively. Darwin or 
Darwin-Riemann ellipsoids are used for the compressible congruent stars without or 
with internal motions. 

For Roche and Roche- Riemann ellipsoids, 

Alyn > ^RL • (4-2) 

This means that for these binary systems secular instability (Roche ellipsoids) or 
dynamical instability (Roche- Riemann ellipsoids) sets in before the Roche limit states 
are reached. 

For Darwin ellipsoids, there is a critical value of the polytropic index as follows: 

r£»2 , (4-3) 

and 

for n < n°, Aiyn > A c , (4-4) 
for n <; n„, A c comes first (no Ad yn ) ■ (4-5) 

This means that for rather stiff equations of state (n < n®.) synchronously rotating 
compressible binary star systems suffer from secular instability and develop internal 
motions. However, for soft equations of state, stars come to contact without suffering 
from any critical phenomena. 

For Darwin-Riemann ellipsoids, the situation is as follows. The critical poly- 
tropic index for the Darwin-Riemann ellipsoids is 



(4-6) 
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and 

for n < n° R , A dyn > A c , (4-7) 

for n ^ ^ct R , A c comes first (no Aiyn) • (4*8) 

Thus for rather stiff binary stars (n < n^) dynamical disruption would occur. 

4.2. Numerical solutions to exact basic equations for binary configurations in New- 
tonian gravity 

Very recently exact equations for binary configurations in Newtonian gravity 
have been numerically solved by us.tP In our treatment, we assume 1) the polytropic 
equation of state, 2) the stationary state in the rotating frame, 3) inviscid fluid and 4) 
irrotational flows in the inertial frame. We do not assume the shape of configurations 
but we solve for the shape of configurations together with other quantities. 

Basic equations are written as follows: 

/ — + (/>+-|W[ 2 - (fi x r)V>F = constant , (4-9) 
Jo 2 



r 



d 6 r , (4-10) 



V6> 

A& = n(H x r - W) • — , (4-11) 



where (9 is defined as 



p = K0 1+n . (4-12) 
Boundary conditions for these equations are as follows: 

0(r s (6,(p)) = for the matter , (4-13) 
n s ■ v s = for the flow , (4-14) 

where r = r s (6,(p), n s and v s are the surface of the star, the normal vector to the 
surface and the velocity vector along the surface, respectively. 

In order to show only the effect of the shape and the structure, we will present 
the results for n = polytropes or incompressible fluids. In Figs.7-^ n results from 
the ellipsoidal approximation EP and those of numerical computations & are shown. 
As seen from these figures, the exact treatment gives larger values of the moment 
of inertia and so the separation at the turning point shifts to larger value compared 
with that of the ellipsoidal approximation. 

For numerical solutions of irrotational Darwin-Riemann compressible configura- 
tions, there is a critical polytropic index as seen in Fig. 9 (eDR means "exact" 
Darwin-Riemann). Therefore, irrotational polytropic binary star systems can be 
classified as follows: 

< r Di? «0.7 , (4-15) 

and 

for n < nf K , A dyn > A OF , (4-16) 
for n ;> rec? R ) ^of comes first (no ^4dyn) • (4 - 17) 
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Fig. 7. The total angular momentum - separation relation for the irrotational Roche-Riemann 
incompressible sequence. The result from the ellipsoidal approximationCP (dashed curve) and 
that of numerical computationsE? (solid curve and dots) are shown. The solid curve and the 
dots are computed by independent numerical methods. 
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Fig. 8. Same as Fig. 7 but for the irrotational Darwin-Riemann incompressible sequence. 




Fig. 9. The total angular momentum - separation relations for several polytropes of irrotational 
Darwin-Riemann sequences. Polytropes with n = 0.0 (dash-dotted), n = 0.5 (dotted), n = 0.7 
(long dashed), n = 1.0 (solid) and n — 1.5 (short dashed) are shown. 
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As discussed before, the critical value n^r « 0.7 is much smaller than the critical 
value 1.2 obtained from the ellipsoidal approximation. This is an impor- 

tant difference because the effective polytropic index for realistic neutron stars is 
estimated 0.5 ~ 1.0. Thus realistic neutron stars may not suffer from dynamical 
instability during evolution due to gravitational wave emission but result in mass 
overflow. 

As for the numerical solutions of irrotational Roche-Riemann compressible con- 
figurations, critical polytropic indices depend on the mass ratio of the stars. 
When n > n^ R , there appear no turning points and so binary star systems would 
result in mass overflow. If n <S n^ R , turning points appear and dynamical disrup- 
tion would occur. In Fig. 10, the critical polytropic index is shown for different mass 
ratios for irrotational Roche-Riemann compressible binary star systems. This result 
is clear contrast to that of the ellipsoidal approximation because solutions from the 
ellipsoidal approximation always give a turning point along the sequences as shown 
in Fig. 11. 

1 
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Fig. 10. The critical polytropic index is plotted against the mass ratio log(Mi/M2) for irrotational 
Roche-Riemann compressible configurations. If the polytropic index is smaller than the critical 
polytropic index, the polytropic star would reach dynamically unstable state and disruption 
would occur. If the polytropic index is larger than the critical value, the mass overflow to the 
other star would begin. 
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§5. Discussion and Summary 

5.1. Discussion 

In this paper we have focused on the equilibrium sequences of binary star systems 
in Newtonian gravity. However, for compact binary star systems such as neutron star 
- neutron star binary systems or neutron star - black hole binary systems, the role 
of general relativity is important. The internal structure and the orbital nature are 
affected by general relativity. Concerning the internal structure, general relativity 
tends to make the matter softer or increase the "effective" polytropic index so that 
the mass will concentrate towards the central region. This implies that mass overflow 
will be more likely to occur. 
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Fig. 11. Same as Fig.7 but for the irrotatioaal Roche-Riemann n = 1 polytropic sequence. T 
results from the ellipsoidal approximation^ (solid curve) and that of numerical computations! 
(dots) are shown. 



As for the orbital motion, from the perturbative analysis of the black hole - 
neutron star binary system, the following relation for the innermost stable circular 
orbit has been obtained: B'0 



d GR = 6M tot + Afi 



where fi is the reduced mass: 



M 1 M 2 



Mi + M 2 

From this radius, we can roughly estimate that 

Mi Mi 
d OF (~ d RL ) > d GK , for — l - Z 0.25, 0.5 £ n £ 1.0 and — ± £ 0.15 

M-2 -ft* 



(5-1) 



(5-2) 



(5-3) 



This means that the mass overflow radius is reached before the general relativistic 
effect begins to work. 

Concerning the critical radius of the dynamical instability, ShibataliiP has shown 
by using the post-Newtonian approximation that 



<^dyn (f° r post — Newtonian case) < dd yri (for Newtonian case) 



(5-4) 



Thus the qualitative results about the critical radii obtained from the Newtonian 
analysis would not be changed even in general relativistic analysis, although the 
values are necessarily different. 

We need to comment on the behavior of the maximum density of the individ- 
ual star along the quasi-stationary evolutionary sequence. In general relativistic 
computations, Wilson, Mathews and MaronettiEJ) have shown that the maximum 
density increases as the system evolves^ However, recently, Bonazzola, Gourgoulhon 
and MarckEj) and Uryu and EriguchitS have shown that the central density of the 
component star of the irrotational binary star systems in general relativity decreases 
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as the binary star system evolves, although at large separations there is some stage 
where the central density increases slightly. 

Our Newtonian computations also show that the maximum density decreases as 
evolution proceeds. This is shown in Fig. 12. For Newtonian polytropes, Taniguchi 
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Fig. 12. The maximum density is plotted against the separation for n = 1 irrotational Darwin- 
Riemann configuration systems. The maximum density is normalized by the maximum density 
of the spherical star. 



and Nakamura@) have analyzed irrotational binary star systems of n 
perturbationally and found that 



- Pn 



45 



Pu 



2{M l /M 2 ye 1 



1 polytropes 



(5-5) 



to the order of 0(1/A e ), where p m ax,i and ao are the maximum density and the 
radius for the initial spherical star, respectively. Here £i equals to ir for n = 1 case. 
Therefore, it is very likely that the maximum mass will not increase significantly. 
In other words, the component star would not suffer from dynamical instability to 
collapse. 

5.2. Summary 

We have reviewed properties of classical and/or approximated equilibrium se- 
quences of Darwin, Darwin-Riemann, Roche and Roche-Riemann problems in New- 
tonian gravity in the context of the final evolution of compact binary star systems 
due to gravitational wave emission. 

We have presented the results of full numerical computations of equilibrium 
sequences of irrotational binary star systems in Newtonian gravity. From the analysis 
of critical points along these equilibrium sequences, we have shown the possibility 
that evolution of compact binary neutron star systems would result in mass overflow 
instead of dynamical disruption for a certain parameter ranges. 
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